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^ ■ Abstract 



Bodeker has recently shown that the high temperature sphaleron rate, which mea- 



> 

| sures baryon number violation in the hot standard model, receives logarithmic correc- 

tions to its leading parametric behavior; T = K'^og{m^ / g 2 T)+0(\)]{g 2 T 2 /m 2 D )a w T A . 
After discussing the physical origin of these corrections, I compute the leading log co- 
efficient numerically; k! = 10.8 ± 0.7. The log is fairly small relative to the 0(1) 
"correction;" so nonlogarithmic contributions dominate at realistic values of the cou- 
pling. 
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It has been known for some time now that baryon number is not a conserved quantity in the 
minimal standard model |L]. It is violated nonperturbatively because of the anomaly, the 
chiral couplings of fermions to SU(2) weak, and the topologically nontrivial vacuum structure 
of SU(2). However, as is characteristic of a nonperturbative process in a weakly coupled 
theory, the rate of violation is so tiny that it is completely irrelevant phenomenologically. 
Certainly, if baryon number is also violated due to high dimension operators descended from 
some GUT, baryon number violating decay rates due to the GUT mechanism will greatly 
exceed the electroweak rate, of order m w exp(— An/ aw) < 10~ 170 GeV even before accounting 
for additional suppression from powers of small CKM matrix elements and a high power of 
(rn piotoa /m w ) . 

However, as first realized by Kuzmin, Rubakov, and Shaposhnikov in 1985, the efficiency 
of standard model baryon number violation is very much higher at finite temperature [@. A 
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perturbative estimate || based on a saddlepoint expansion about Klinkhamer and Manton's 
sphaleron indicated that the rate is more than enough to erase any relic abundance of 
baryon number left over from the GUT scale, unless the combination B — L, baryon minus 
lepton number, is nonzero; this quantity is preserved by electroweak physics. It also may 
be that the baryon number abundance in the universe was produced by electroweak physics, 
which motivates the more careful study of electroweak baryon number violation at finite 
temperature. 

Recently our understanding of thermal baryon number violation has improved, though it 
is not complete. We understand that the efficiency of baryon number violation can in most 
relevant circumstances be related by a fluctuation dissipation argument to the diffusion 
constant for Chern-Simons number [||, |6|, [7]. This diffusion constant, called the sphaleron 
rate, is defined as 

rsBm h«MiM)n, (1 ) 

V—oot-XJO Vt 

where the expectation values refer to a trace over the equilibrium thermal density matrix. 
The quantity T is the topological susceptibility of the electroweak sector at finite tempera- 
ture, in Minkowski time. 

We know that, when the electroweak phase transition is first order, the value of T jumps 
discontinuously between the phases; it is quite small in the broken phase and much larger 
in the symmetric phase. An old power counting argument says that the symmetric phase 
rate should be of order T ~ <y w T 4 with an order unity coefficient. This argument relied, 
correctly, on the natural nonperturbative length scale in the hot plasma being l/(awT). One 
then assumes that the natural time scale is the same; on dimensional grounds the spacetime 
rate of Nqs diffusion must then be of order (awT) 3 awT. However, Arnold, Son, and Yaffe 
have shown that at leading parametric order, the natural time scale is not l/(awT), but 
l/a w T [pfl, up to possible logarithmic corrections, which the authors did not consider. More 
recently, Bodeker has demonstrated that logarithmic corrections to their argument do occur 

This says nothing about the numerical value of the sphaleron rate; it could be paramet- 
rically a w T A but numerically irrelevantly small. Ambj0rn and Krasnitz presented numerical 
evidence that it was large, by considering classical, thermal Yang-Mills theory on the lattice 
|T0| . Their definition of iVcs was not topological and could therefore suffer from potentially 
severe lattice artifacts, but Turok and I studied the same system with a topological defini- 
tion of iVcs and verified that T is substantial [|ll| . The sphaleron rate we found, expressed 
in physical units, was lattice spacing dependent, which turns out to be a prediction of the 
arguments of Arnold, Son, and Yaffe; the extra power of aw arises from the interaction 
between the infrared fields and ultraviolet excitations, as I will discuss more below, and on 
the lattice it becomes one power of the lattice spacing a. 

In fact the lattice spacing dependence of T for pure classical lattice Yang-Mills theory 
only fits r oc a if there are substantial corrections to scaling, which has led some to call 
into question whether the Arnold, Son, and Yaffe's analysis is correct. In this paper I will 
assume that it is, as seems justified on theoretical grounds and numerical evidence from 
classical Yang-Mills theory "enhanced" with added degrees of freedom which reproduce the 
hard thermal loop effects [0 . I will return to the large corrections to scaling in the classical 



lattice theory in subsection [5]3. 
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Because the interactions between infrared and ultraviolet excitations are important to 
setting the sphaleron rate, it is somewhat difficult to actually extract T at leading order in 
aw for the continuum quantum theory, at the physical value of aw or even in the parametric 
small aw limit. Hu and Miiller proposed a technique based on including the UV physics, lost 
to the lattice regulation, by introducing new degrees of freedom which influence the IR fields 
in the same way ||12|| . They implemented and applied the technique jointly with me [|l3j. As I 



will discuss later, this technique still suffers from some poorly controlled systematics, which 
are related to the logarithmic corrections discovered by Bodeker; in fact, beyond the leading 
log the infrared physics the technique will simulate is not rotationally invariant. 

A good first step to answering the remaining questions about the sphaleron rate is to 
determine it at leading logarithmic order in aw- Bodeker has demonstrated that this can be 
done within an effective theory which is completely UV well behaved; in fact it is nothing but 
the Langevin equation for 3-D Yang-Mills theory ||. The leading log behavior is probably 
not very useful by itself, for estimating T at the physical value of aw- Leading log expansions 
often miss large constant corrections; we know for instance that the 0(g 2 T) contribution to 
the Debye mass has a much larger constant contribution than the leading log might suggest 



|14] , [15]. However, it is still useful to know the leading log; for instance, its size is related 
to the severity of the systematic problems with the method of Hu and Miiller, and it might 
in principle be useful for extrapolating lattice results which correspond to an inappropriate 
value of a w back to the correct value. 

The purpose of this paper is to determine the coefficient of the leading log behavior of 
T; namely, to find k' defined through 

r = «' f log ^ + 0(1) J ^a 5 w T A + (higher order) . (2) 



g 2 T J m 



(My logs are always natural logs.) Neglecting the 0(1) means that log(l/g) is treated as much 
larger than any order unity constant, a rather extreme interpretation of the perturbative 
expansion. Probably this expansion is completely unjustified, but as I said there are still 
important things to be learned from making it. Until Section [| I will not worry about 
whether the expansion in log(l/g) ^> 1 is justified; the goal is simply to determine k' . I will 
also work in Yang-Mills theory, which is appropriate at leading order only for temperatures 
well above the equilibrium temperature. I will mention how to include the Higgs field in 
Section 

A summary of the paper is as follows. Section || will discuss in an intuitive, physical, but 
nontechnical level why the rate has the parametric form I show; where the extra a comes from 
and particularly why there is a log. The section provides two apparently different arguments, 
one in terms of conductivities and scattering processes for hard particles and one in terms 
of hard thermal loops (HTL's) and Wilson lines; the two are of course equivalent. Some 
salient details about the Wilson line are in Section |3|, which is more technical. This section 
also shows how the numerical model of Hu, Miiller and Moore fails beyond leading log order. 
Section |] studies Bodeker's effective theory for extracting the leading log, numerically. Since 
the Langevin dynamics are simple but numerically costly, the emphasis is on controlling 
systematic errors. Many of the details appear either in previous papers or in the appendix. 
For the reader's ease I present the answer now: k! = 10.8 ± 0.7. The dominant error here is 
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statistical; systematic errors are completely under control. In Section |5| I discuss the meaning 
of the result, and try to estimate what the sphaleron rate is for the realistic values of m 2 D 
and g 2 by using this result to extrapolate the results of Hu, Miiller, and Moore to the right 
value of the log. My estimate is T ~ 20a^/T 4 for m 2 D = (ll/6)g 2 T 2 and g 2 = 0.4; however 
there are uncontrolled systematics which may be as large as 30%. I also discuss corrections 
to the approximation that g <C 1; the largest of these will be of order 10%. The conclusion 
concludes. There is also a technical appendix which discusses the match between lattice and 
continuum Langevin time scales; the match between continuum and Langevin time scales 
is computed at 0(a) and the match between lattice Langevin and heat bath time scales is 
determined by a measurement. 



2 The physics behind a\y log(mz)/g 2 T) 
2.1 note on the classical approximation 

To a very good approximation the behavior of infrared fields in thermal Yang-Mills or Yang- 
Mills Higgs theory at weak coupling is that of classical fields. The "old" parametric estimate, 
T oc a^T 4 , relies on the fact that the only length scale available in classical Yang-Mills theory 
is l/g 2 T. It is also assumed that the only time scale is the same, in which case any essentially 
infrared spacetime rate must go as g 8 T 4 . 

This argument relies on a decoupling between the infrared and ultraviolet fields, since the 
ultraviolet fields (by which I mean fields of wave number k ~ T) do not behave classically. It 
is known that this decoupling is very accurate for thermodynamic variables, except that the 
Aq component of the gauge field receives a Debye mass. In fact the validity of the decoupling 
is equivalent to the quality of the dimensional reduction approximation [|16| , which has been 
discussed extensively ^T7\, [THJ. 

However, the decoupling does not extend to dynamics; the generalization of the Debye 
mass to unequal times are the hard thermal loops, which significantly affect the infrared 
dynamics. The hard thermal loops are precisely that set of diagrams which are linearly 
divergent within the classical theory. This linear divergence is cut off at the ultraviolet scale 
k ~ tcT where the theory ceases to behave classically; so the size of these effects depends 
essentially on the way the IR classical theory is regulated.^] Thus it is only correct to say that 
the IR fields behave classically if we mean that they behave like the IR fields of a classical 
theory regulated in some way which correctly reproduces the hard thermal loops. In nature 
that regulator is quantum mechanics, but we might be able to find some other appropriate 
regulator in a numerical setting. 



2.2 Argument in terms of conductivity and scatterings 

Here I will give an argument for the a^y log(l/g)T 4 law based on Lenz's Law and the con- 



ductivity of the plasma. The argument has very recently been made quantitative fl9j| , but I 
will present it at the qualitative, intuitive level. 



*I should mention parenthetically that the classical theory also contains quadratic and cubic divergences 
in the energy density, but these do not affect the IR dynamics responsible for baryon number violation. 
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To see how hard thermal loops influence the sphaleron rate, I first make the point that the 
sphaleron rate is set by the evolution of very soft infrared fields, where by very soft I mean 
fields with wave number k ~ g 2 T. Parametrically shorter wave lengths do not contribute 
appreciably because the probability for nonperturbative physics to occur at such scales is 
exponentially suppressed. According to standard sphaleron type arguments, the contribution 
from the scale k ~ g 2 ~ s T is suppressed by of order exp(—k/g 2 T) = exp(—g~ s ). Even the 
scale k ~ g 2 T log(l/g) gives a contribution suppressed by a power of g. 

The second point is that it is physics in the transverse sector which matters, and in 
particular, diffusion of Acs requires the evolution of magnetic fields. To see this, first go 
back to the definition of Acs,Q 

Ncs(h) - Ncs(h) = f 2 dt [ d 3 x^E«B»(x) . (3) 

Now the B field is always transverse, meaning that D ■ B = 0, by the Bianchi identity; and 
so only the transverse part of the E field contributes. The Bianchi identity also states that 
D x E = —[Dt, B\. Since the relevant part of the electric field is transverse, it will in general 
have nonzero covariant curl. For E to be nonzero and to remain the same sign for long 
enough to give a nontrivial contribution to / EfB^dt, there must then be time evolution of 
infrared magnetic fields. 

At this point it is useful to recall how infrared magnetic fields evolve in the abelian 
theory, on wave lengths longer than the Debye screening length. The answer is familiar 
plasma physics; the plasma is very conducting, and a conducting medium resists changes 
in magnetic fields by Lenz's Law. In the limit of infinite conductivity the magnetic fields 
are perfectly frozen; for finite conductivity the time scale for their evolution scales with 
conductivity. A magnetic field of wave number k, with well more than its mean thermal 
energy density but much less energy density than is contained in the bulk plasma, decays 
according to 

d A ■ 1 

^ = ^ = -7Ik^) {k% '- k ' k ' )A " (4) 

in the parametric limit that the decay time is well longer than 1/k, which is satisfied at 
all length scales parametrically longer than the Debye screening length. The characteristic 
decay time of a magnetic field in the plasma is then r = a(k,uj <C k)/k 2 . Of course I have 
only written the dissipative part of the magnetic field evolution equation; there must also 
be a noise term which is uniquely specified by the requirement that the thermodynamics of 
the IR magnetic fields are correct. 

Note that the conductivity is wave number dependent. It has a good infrared limit which 
is achieved for length scales larger than the mean scattering length /f ree of a current carried 
by hard particles. In the abelian theory a particle's charge is preserved when it undergoes a 
scattering, so this length scale is the mean length for large angle scattering, parametrically 
Zfree ~ 1 / a 2 T times logarithmic corrections. If every particle had the same free path and a 
scattering perfectly randomized its momentum, the conductivity on scales longer than Zf rcc 

2 Although I set the speed of light c = 1, I typically write expressions noncovariantly with a positive 
space metric, which is convenient in the finite temperature context because the thermal bath establishes a 
preferred frame. 
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would be 

a = , k « ■ (5) 

When scattering processes are more complicated this formula defines an effective value of 
/free- For length scales well between the Debye screening length and Zf ree , where scatterings of 
the charge carriers can be neglected, the conductivity is related to the Debye length through 

a = —f- , m D > k > 1/Zfree • 6 

4k 

(In both expressions the Debye length mo, which is O(gT), is just keeping track of the 
number density, charge, and (l/E) of the particles. The derivations of each expression assume 
ultrarelativistic dispersion relations for the charge carriers.) This expression follows from the 
form of the transverse self-energy and the fact that the conductivity we are discussing in this 
case is just a special case of the HTL self-energy, a(k,uj <C k) = ImU T (k,u})/u. This is the 
connection between this "conductivity" picture and the hard thermal loops. 

The behavior of the electroweak gauge fields for scales parametrically between the non- 
perturbative scale k = g 2 T and the Debye scale k = gT is the same as in the abelian theory 
at leading parametric order.0 Since the mean free path of a hard excitation to undergo any 
scattering is ~ 1 / (g 2 T log(l/g)) |21f , the conductivity for k = g 2 ~ 5 T, with < 5 < 1, is 
a = nm 2 D /4k ~ g s T, and the time constant associated with the decay of a magnetic field 
is t ~ l/(g 4 ~ 3<5 T). Although the scale k ~ g 2 T does not fit within the range of validity of 
this argument it cannot be that the decay rate for a magnetic field with k ~ g 2 T differs 
from the 5 — > limit by any nonzero power of g. Hence the relevant time scale for the 
dynamics of nonperturbative infrared magnetic fields in weakly coupled, hot Yang-Mills the- 
ory is 0((g 4 T) -1 ), up to corrections at most logarithmic in g. This is a paraphrase of the 
argument of Arnold, Son, and Yaffe ||, who however neglected the possibility of logarithmic 
corrections. 

Logarithmic corrections do in fact occur. While in the abelian theory the electrical 
conductivity only reaches a long wave length limit at a scale set by a hard particle's free 
path for large angle scattering, in the nonabelian theory that limit is set by the free path for 
any scattering. The reason is that when a colored particle undergoes a scattering, however 
small the transfer momentum, its color is changed. This degrades the color current even if it 
does not degrade the momentum carried by the particle substantially. The total rate for any 
scattering to occur is twice the damping rate, which has been computed for hard particles 
at leading log by Pisarski. The damping rate for an adjoint charged particle of any spin is 



21 



_ Ng 2 T 
7 4tt 



log ^ + 0(1) 
9 T 



(7) 



and for a fundamental representation particle it is the same with N — > (N 2 — 1)/2N. Note 
the log, which arises from an integral over exchange momenta, running from the gT to the 

3 Some years ago Ambj0rn and Olesen argued that nonabelian fields obey an anti-Lenz's Law ^o), ap- 
parently in contradiction to the argument presented here. Their work refers to the nonabelian interactions 
between W and Z fields in the presence of strong (clcctromagnetism) magnetic fields at zero temperature. 
It may have some bearing on the mutual interactions of the fields at the g 2 T scale, but the response of the 
harder modes to the very soft fields is at leading order the Lenz law type behavior seen in the abelian theory. 
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g 2 T scale; also note that the result is independent of the particle's momentum, provided it is 
harder than the gT scale. A collision largely but incompletely randomizes a particle's charge, 
and so on scales longer than ~ l/(g 2 T log(l/g)) the electrical conductivity of the plasma 
is ~ m 2 D /2> r ), up to nonlogarithmic corrections. To find the numerical constant one must 
determine how thoroughly a scattering randomizes a particle's charge, which depends on the 
representation of the particle. Also, the particle's charge is not destroyed, just transferred to 
another particle; one must check whether this induces any important currents. (It turns out 
not to, because a particle is as likely to scatter from a charge carrier moving in one direction 
as in the exact opposite direction.) The calculation is quite nontrivial but it has been done 
recently by Arnold, Son, and Yaffe [19|], who show that, at leading log, the conductivity is 
simply m 2 D /3 , ~f, with 7 given in Eq. (|7[). This is independent of the group representation of 
the particles carrying the current. 

Performing an extreme parametric expansion, log(l/g) ^> 1, the scale set by particle mean 
free paths and the nonperturbative scale are well separated, and the fields with k ~ g 2 T see 
a k independent conductivity. Hence the relevant infrared dynamics for transverse modes is, 
at leading logarithmic order, 

3Ng 2 T\og(l/g) ^ 
Po. M = 4^ D o F oi + noise > (8) 

where the transverse part of the noise is fixed by the requirement that the thermodynamics 
come out right. The longitudinal part of the noise generates time dependent gauge rotations 
of the A fields, which are irrelevant to Chern-Simons number; so we may choose the amplitude 
of the longitudinal part of the noise to be whatever we want. It is most convenient to choose 
it to be of the same magnitude as the transverse part, in which case we reproduce the effective 
theory of Bodeker, which is also the Langevin equation for 3-D Yang-Mills theory. 

While this derivation has presented the ideas in an intuitive way it is scarcely rigorous, 
so I will now approach the problem a little more formally by looking at the hard thermal 
loop effective theory. 



2.3 HTL effective theory and the Wilson line 

Now I will go through the argument for Bodeker's effective theory also from the point of view 
of the effective HTL theory for the infrared modes. The idea of the HTL effective infrared 
theory is that one can construct an effective theory for the modes with k -C T, valid at 
leading order in g, by integrating out all degrees of freedom with k ~ T. The procedure for 
separating the degrees of freedom and regulating the effective theory, cut off at some scale 
between gT and T, is left unspecified and does not affect the result at leading order. 

The resulting effective theory is a classical theory for the remaining modes, but with a 
nonlocal "HTL" effective action correction. The nonlocal effective action was first derived 
in [22 1 . In the current context, since in the long term we have in mind a position space 
regulation of the IR theory, it is most convenient to write the effective action in real space, 
as was first done by Huet and Son p3f : 

2 

-[D ,E t r(x,t) + [D J ,F Jt } a (x,t) = %(x,t) + ^ Jd 3 y^U ab ((x,t),(x + y,t-y))x 

E b Ax + y,t-y). (9) 
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Figure 1: "It's quite simple, really . . . ." A scorecard of the scales involved in the problem 
and the approximations which are valid in each. None of these scales are distinct if we do 
not take g <C 1. 



Here U ab ({x, t),(x + y,t — y)) is the adjoint parallel transporter along the straight, lightlike 
path between the points (x,t) and (x + y,t — y). Note that the electric field on the right 
hand side is at the retarded time t — y. The noise £° is Gaussian with a nonvanishing two 
point function for lightlike separated points, 

t)$(x + y, 0> = ^K\t - t'\ - \y\) yj ^U ab ((x, t),(x + y, t')) . (10) 

For k <C gT, the effective action simplifies somewhat. The coefficient of the term on the 
right in Eq. @, which contains one time derivative (since E{ = [Dq,Ai]), is large enough 
that the time evolution is determined by this term and we can neglect the [Do, E] term 
on the left. Further, the time scale associated with the fields' evolution is longer than the 
time retardation appearing in the nonlocal HTL action, and we are justified to neglect that 
retardation. The result is 

[Dj, Frfix) = &{x) + ^ / d 3 y^U ab (x, x + y)E%x + y) , (11) 

and the noise correlator is simplified by dropping the retardation there as well (though to get 
its magnitude right we must remember that there are two contributions from 5(\t — t'\ — \y\)). 
If in addition k 3> g 2 T\og{l/g) then it is possible to choose a gauge such that the parallel 
transporter is close to the identity and at leading order it can be ignored. In this case we 
can recover Eq. (^) by expanding [Dj, Fji] to leading order in A and Fourier transforming. 

Huet and Son argue that, to model the 0(g 2 T) modes alone, we can use Eq. flTlD , 
interpreted as an expression for fields at the g 2 T scale alone with all shorter wave length 
scales integrated out. The reasoning is that the higher modes only appear in the equation 
in [Dj,Fji\. But [Dj,Fjj\ is a thermodynamical quantity depending only on the transverse 
fields, and as already stated, the UV causes perturbatively small corrections to this sector, 
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which can be ignored at leading order. The power counting arguments are laid out explicitly 
in M. 



However, this misses one key issue. Can we neglect the influence of the modes with 
g 2 T <C k <C gT on the adjoint parallel transporter U ab (x,x + y) when y ~ l/(g 2 T)? The 
answer is, no. When evaluating an adjoint parallel transporter for a path of length ~ 
l/(g 2 T), essential contributions arise from all scales intermediate between gT and g 2 T. This 
problem has recently been addressed by Arnold and Yaffe, in the context of studying 0(g 2 T) 
corrections to the Debye screening length J25|. They show that, in SU(iV) pure gauge theory, 
any two point correlator of equal time adjoint fields at separation y, connected by a straight 
adjoint Wilson line, falls off at least as rapidly as 

exp(-y/A), A- = ^(bg^ + K), (12) 

where UV regulator and K is a constant of order unity. The constant is evaluated 



for SU(2) in | 15(| ; for the regulation appropriate to determining the Debye mass, m reg = mr, 
and K = 6.7±.3. Nonlogarithmic corrections are large in this case; I will return to this point 
later. The value of the constant K includes nonperturbative physics at the g 2 T scale, but 
the logarithm arises perturbatively from scales intermediate between g 2 T and m reg . Thus, 
to study physics on the length scale l/g 2 T, we must include the influence of scales between 
gT and g 2 T on the parallel transporter in Eq. (|l|). 

Note the sneaky way the modes with g 2 T <C k <C gT have entered the dynamics of 
the g 2 T modes. Their direct influence on the interactions between the g 2 T modes is indeed 
small, as Son showed ||24j| ; but they change the way that the g 2 T modes interact with the 
hard modes, in a way which turns out to be important. Bodeker emphasizes this viewpoint 
in his original derivation of the effective theory for the leading log ||. 

Using the correction to the Wilson line from the intermediate momentum modes, I can 
now establish Bodeker's effective theory for the k ~ g 2 T modes. The modes with k ^> g 2 T 
change quickly compared to the g 2 T scale modes, so the g 2 T fields see an average over the 
more UV scales. In particular the parallel transporter relevant for the g 2 T modes' evolution 
is the average of the parallel transporter over realizations of the g 2 T <C k -C gT fields. 
As I demonstrate in subsection |3.1| , averaging over realizations of k ^> g 2 T modes leads to 
exponential damping of the parallel transporter, for y of order l/(g 2 T\og(rriD/g 2 T)). In the 
Coulomb gauge,(] 

Coulomb x + y) * 5 ab exp(-y/A) , (13) 

with A the same as in Eq. flTSp. If we are permitted to expand in log(m£>/g 2 T) ^> 1, then the 
integrand on the right hand side of Eq. (|11]) has already fallen away before y comes on order 
l/(g 2 T); therefore the approximation which gives Eq. (|3~3|) is valid, up to log corrections, 
throughout the range of y which dominates the contribution to the integral. The effective 

4 The use of Coulomb gauge becomes problematic when considering length scales I > l/g 2 T, and for 
considering any unequal time correlator if the total volume of space considered is V 3> (g 2 T)~ 3 ; however 
this is not relevant because technically we are only applying Coulomb gauge to modes with k 3> g 2 T in order 
to integrate them out and establish an effective theory. The gauge fixing of the 1R effective theory, ie of the 
problematic k ~ g 2 T modes, has not been specified. Also, at leading log the same results would be obtained 
in Landau gauge. 
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theory for the g 2 T modes (fixing the gauge freedom on scales more UV than g 2 T to Coulomb 
gauge) is therefore 

+ / d 3 ye-W^ E «( x + y) = [Dj, F^x) . (14) 

y 

Since the integral is dominated by y ~ l/(g 2 T\og(m D /g 2 T)), and for the IR fields of 
interest E varies only on the l/g 2 T scale, it is permissible at leading log order to pull the E 
field out of the integral. The integral is then quite simple; performing it gives 

[Dj, Fji] = -fE? + £ , (£?(x, t)$(y, t')) = -^5{x - y)5(t - f , (15) 

where the form for the noise correlator also follows from the approximation for the parallel 
transporter. (Alternately, one can always recover the form of £ by insisting that the ther- 
modynamics come out correctly. ) This is Bodeker's effective theory, though it remains to 
establish that m rcg should be m D . I discuss this more in the next section. 

3 Wilson line, more carefully 

Now it is time to look more carefully at the Wilson line appearing in the last section, first to 
verify the claims there, second to show the connection to the argument involving scatterings 
of the hard particles, and third because it is relevant to the analysis of the results of Hu, 
Muller, and Moore. 

3.1 Exponential decay of the Wilson line 

What we want to know about is the Wilson line between lightlike separated points a spatial 
distance I apart, with / 3> l/m D but less than l/(g 2 T) by at least a logarithmic factor. 
Actually, we want to know the average of the Wilson line over realizations of the modes with 
k > g 2 T\og(l/g), since these fields fluctuate faster than the 0(g 2 T) fields, which therefore 
see the average over realizations, up to corrections subdominant in l/log(l/g). 
The Wilson line is given by 

U = Pexp[ ig(A + A z )(z,t = z)dz , (16) 
Jo 

with A in the adjoint representation. I will only consider the transverse contributions here, 
the longitudinal ones are subdominant. I use Coulomb gauge, in which the Aq contribution 
arises entirely from the longitudinal modes and the A z only arises from the transverse modes. 
Further, I will work at leading parametric order, by which I mean that higher point correlators 
and vertex insertions are ignored, and combinations of A fields are evaluated assuming A 
is Gaussian by applying Wick's theorem. Of course I will include the hard thermal loop 
corrections to the A field propagators. These approximations are justified at leading log 
down to k ~ g 2 T\og(l/g), which is all we need. (If I were interested in lengths I ~ ^-/g 2 T 
rather than logarithmically shorter, most of the approximations I make would break down 
completely.) 
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All terms with odd powers of A vanish on averaging over realizations, while even terms 
look like (applying Wick's theorem between lines 1 and 2) 

oo / \n n 2n „ 

( u ) = E Sfrr / dZl ■ ■ ■ dz m(Mzi) ■ ■ ■ M**)) (17) 

n=0 \ ZU )- J 

00 1 ( -a 2 r l V 
= El "f dz.dz^A^M^))) 
^0 n] V 2 Jo J 

1 2 rl 



exp (-^~ J dz 1 dz 2 (A z (z 1 )A z (z 2 )} 



so the average over realizations is the identity times the exponential of the two point con- 
tribution. (Note that for any given realization the Wilson line has unit modulus. But the 
average over realizations does not, its modulus falls exponentially with distance.) 
To evaluate this we need the two point correlator, 

(Mx,t)My,t>)) = J -0^I p (*,a,) ^ - e ^) e -Mv^ . (18) 

Here p(k, u) is the spectral density, which is the magnitude of the discontinuity in the 
propagator 1/(uj 2 — k 2 — IIy(A;, u)) across the real to axis on analytic continuation from 
Euclidean (imaginary) uj, 

p(k, uj) = 21m Uu + it) 2 - k 2 - U T (k, uj + ie)) ^ , (19) 



and T/uj is the classical approximation for l+n(u), with n(u) the Bose distribution function. 

So far I have suppressed group indices, but when they and the integrals over z are 
evaluated we get 

Tr (U) Ng 2 r d 3 k du Tp(k, uj) k 2 x + k 2 y 4 sin 2 (A; 2 - u)l/2 
° g ^r7l~~ T 1 J2nf2^ u k 2 (k z - uj) 2 ' ( ' 

The equal time correlator of the A field goes as T/k 2 plus subleading corrections, for all 
k > g 2 T; so 

du Tp(k,uj) _ T 

2tt uj ~ k 2 ' [ ] 

What matters now is where p is concentrated. In the regime k -C mo, almost all the 
contribution to Eq. ([H]) is from \u\ k, see [^T|. This is just the statement that these 
modes evolve on time scales slower than 1/k, as I have already discussed. Hence, in evaluating 
the low k contribution to Eq. ( p0|) I can set k z — uj ~ k z , and then perform the integral over 
uj, giving 

Ng 2 T fW<m D d 3 k 1 kl + k 2 4sin 2 /U/2 

(2^P~P k 2 ' [22) 



The k z integral is completely well behaved; in fact the large / limit of 4 sin 2 (k z l/ '2) /k 2 is 
2irl5(k z ). This represents the fact that only modes with k z < l/l have A of the same phase 
all along the Wilson line; other modes' contributions destructively interfere in the integral 
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along the Wilson line. The integral over the other two directions is dominated by a logarithm 
arising from scales intermediate between k± ~ l/l and k± = mp, where the approximation 
for p breaks down. Performing the integral over k± first gives 

Ng 2 lT ( m D \ f sm 2 (k z l/2) MU Ng 2 lT ( m D \ 

In the opposite limit, k 3> mo, the excitations obey normal vacuum ultrarelativistic 
dispersion relations to good approximation, so the spectral density is approximately 

p(k^m D ,u)~-(8(k-u)-5(k + oj)). (24) 
k 

The large k contribution is then 

Ng 2 r d 3 k k 2 x + k 2 y T 4sin 2 ((|A;| -k z )l/2 



(2tt) 3 k 2 k 2 (\k\ - k z ) 2 



(25) 



Now the term with sin 2 in it is forcing \k\ = k z ; only modes propagating along the direction 
of the Wilson line keep in phase, others destructively interfere. But the polarization vector 
of such a mode is close to orthogonal to the Wilson line; because k ^> l/l, for any k for 
which \k\ — k z < l/l, the factor k\/k 2 will be near zero. Thus, the only modes which avoid 
destructive phase interference are polarized in the wrong direction to contribute significantly. 
Continuing to carry out the integral by defining x = k z /k, we get 

N 9 2 T f 1 , f „, 2 Asm 2 (lk(l-x)/2) Ng 2 T 
-~fr dx dk(l-x 2 ) , V >' > ~ log (lm D )+0 1 , 26 

Z7T Z J-l Jm D (1 — Xy 2lT z m£) 

which is 0(g). The hard modes do not contribute at leading order to the Wilson line along 
a lightlike path. 

I have not treated the modes with k ~ mo, which are more complicated because this is 
where p does not fit into either limiting category. But they turn out to give a result smaller 
by a logarithm than that from the modes with fc<mj). The final result is that, on averaging 
over realizations of the modes with k greater than g 2 T by at least a logarithm, the Wilson 
line in Coulomb gauge is 

(U) = 1 expH/A) , A" 1 = ^ (log p| + 0(1)) . (27) 
3.2 relation to scattering 



Now look at Eq. ( pO|) again. Remember that the Wilson line is representing the trajectory 
of a hard particle, with p much greater than any momentum scale which gives a leading 
order contribution to the integral. We want to interpret Eq. ( P0|) as I times the rate for the 
particle to undergo a scattering involving the transfer of a soft field, times a (representation 
dependent) group theory factor which tells how thoroughly the scattering randomizes the 
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(a) (b) (c) 

Figure 2: (a) Emission of a gluon below the light cone; (b) Self-energy diagram which, cut, 
describes the emission; (c) The cut must split an HTL insertion, so the physical process is 
scattering with a small exchange momentum. 



particle charge. To see the relation, take the large / limit. Then 4 sm 2 ((k z — u)l/2) / {k z —oj) 2 = 
2irl5(k z — to), and the decay rate per unit length of the Wilson line is 

Ng 2 r d 3 k du Tp(k,uj) k 2 x + k 2 

— Jj2^2;— jfci-^to-")- (28) 

Recall that the hard particle starts out with p z ^> \k\, p x — p y — 0. If it emitted a particle 
of wave number k, its momentum would change to p — k and its energy would change by 
—k -p = —k z plus a correction of order k 2 /\p\, which is negligible by assumption. The delta 
function is just the energy conserving delta function appearing in the expression for the rate 
of the process shown in Figure (a), the emission of a soft gluon by an adjoint charged, hard 
mode. 



The rate for this process has been considered by Pisarski [21]; his expression (4.4) does 
not look quite the same as Eq. (|25p for the following reasons. First, his definition of the 
spectral density differs from the one used here by a factor of 2ir. Second, he has allowed 
the hard particle to be a little off shell, which is important in the part of the integral with 
k ~ g 2 T log(l/g) but not for higher k. (Note however that my large I approximation breaks 
down in exactly this regime.) For momenta k 3> g 2 T log(l/g) we may integrate over uo t in 
his expression and obtain Eq. fl2"8|). The difference caused by k ~ g 2 T\og{l/g) does not 
change the leading log, but would be important in investigating corrections to the leading 
log. Also note that Pisarski considers the contribution of longitudinal momenta and finds 
that they do not give a log, only a constant times Ng 2 T. 

The process considered involves the emission of a soft mode with \k\ > \u\. The reason 
that such emission is possible at all is because of the hard thermal loop correction to the 
gauge propagator. The rate is the imaginary part obtained by cutting the one loop self-energy 
diagram, Fig. 0(6). The cut gives a nonzero result when it goes through a HTL self-energy 
insertion, which can be interpreted as diagram (c) in the figure. Hence the Wilson loop 
calculation is proportional to the rate for scatterings of one hard particle off all other hard 
particles by exchange of soft intermediaries. 

Pisarski calculates the total rate for a particle to undergo a collision, and because of 
the conventional definition of the damping rate, his damping rate is actually half of the 
scattering rate. Also, his result does depend on the representation of the particle undergoing 
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the scattering. The rate the Wilson line calculation determines is actually the rate of color 
randomization, not of collisions; there is a representation dependent correction between 
the two, which depends purely on color factors at the vertex where the particle of interest 
interacts with the soft background field. The total collision rate of a particle is proportional 
to the group factor 

<~r\. r pa'Y [ a 

(») 

where T a is in whatever representation the particle is in. The original disturbance of the 
particle distribution from equilibrium is caused by an electric field, which is an adjoint 
object; the disturbance to the single particle density matrix is proportional to E b T b . The 
color randomization per collision is 

Trl T a T b T a . . 

(30) 



where the second term tells how much the color after the collision is aligned with the color 
before. Multiplying by the total collision rate gives a color randomization rate oc N/2, 
independent of representation. 

3.3 subleading corrections in the method of Hu and Miiller 

Now I will discuss the relation between the quantum theory and the technique proposed by 
Hu and Miiller, refined and implemented jointly with myself, concentrating on whether the 
behavior is the same at next to leading logarithmic order. 

First I should explain why the technique is necessary. Traditionally people have tried to 
determine the sphaleron rate by studying classical Yang-Mills (or Yang-Mills Higgs) theory 
regulated on a spatial lattice. However, if we just study the classical theory on the lattice, 
using the lattice spacing as the UV cutoff, the hard lattice modes generate HTL effects which 
do not look like Eq. (|9|) and in fact are not rotationally invariant [p6| . Even in the leading 
log approximation this is a problem, because it means that, where Eq. flT5|) has E{, we will 



get E rescaled by a rotationally non-invariant factor, determined in |27 . 

A proposal by Arnold to fix this problem, staying within lattice classical theory, by 
making the dispersion relations turn up very steeply [p7[| , does not work because the hard 
modes are then Landau dampedP]. The only alternatives I am aware of involve adding new 
degrees of freedom which influence the IR classical fields in a way equivalent to correct hard 
thermal loops. Two such proposals exist in the literature. One is due to Bodeker, McLerran, 
and Smilga p6| , more recently discussed by Iancu p8| . I will not discuss it since no one has 



yet specified a complete discrete numerical implementation. The other idea was proposed by 
Hu and Miiller [|12[]; the details of the implementation were worked out and applied jointly 
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We also simulated the classical IR physics by studying nonperturbatively the classical 
system regulated on a spatial lattice, thereby treating the left hand side of Eq. (P) fully 



The problem with Arnold's proposal is that it gives a Lorentz non-invariant hard mode dispersion relation 
under which 1 — > 2 and 2^1 processes are kinematically allowed. They are efficient, so the hard excitations 
have a mean free path for hard scatterings ~ l/g 2 T. I have Arnold's agreement on this point. 
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Figure 3: A simple 2-D example of how a Wilson line (diagonal) on a lattice (dashed lines) 
is replaced by the sequence of links which stay closest to it (solid, jagged line). 



nonperturbatively. To include the HTL effective action, the right hand side of the expression, 
we added to the classical lattice system a large number of adjoint charged classical particles. 
They take coordinate positions in the continuous space in which the lattice fields sit, obey 
ultrarelativistic dispersion relations, and interact with the lattice fields when they cross the 
dual planes to lattice links. There are two parts to the interaction with the lattice fields. 

First, the particles "kick" lattice electric fields, and their momenta receive a similar 
"kick". The kicks and the approximately random distribution of the charges performs the 
noise, and a correlation between past fields and the "kick" the gauge field receives, arising 
from a change in the particle trajectory from the "kick" it received, accounts for the nonlocal 
term. The size of each kick is proportional to a charge Q which is made small so the particles 
individually interact weakly. 

Second, the particles' adjoint charges are parallel transported by the gauge field connec- 
tion. Since the gauge fields are only defined on the lattice links, the Wilson line used for 
the parallel transportation of a particle is the sequence of lattice links which maintains the 
shortest distance to the actual path of the particle, as illustrated in Figure [|. The exact 
rule used to choose the sequence of links is that a link is used if the Wilson line penetrates 
the face of the dual lattice corresponding to that link. Parallel transporting the particle's 
charge ensures that we are including the parallel transporter in the right hand side of Eq. 
(§); it is also absolutely necessary to make the update rule manifestly gauge covariant. The 
full update is described in excruciating detail in [13| . 

The Debye mass depends on the density (n) and charge Q of the added particles as 
m 2 D ex Q 2 (n), so by changing the number of particles one can tune the HTL strength to be 
whatever is needed while keeping Q small. If there were only IR fields, meaning fields with 
k <C 1/a (a the lattice spacing) then the behavior would correctly reproduce Eq. (H), at 
least in the limit Q — > 0, (n) — > oo with m 2 D fixed. This is also discussed in |L3 . 



The lattice theory involves two scales, a and mr>\ to get the leading parametric behavior 
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we must seek the limits a — > 0, nip — > oo (if we think of the length scale l/(g 2 T) as 
remaining fixed). There are two ways we could go about doing this. One corresponds to 
using a as a cutoff between the gT and the T scales, meaning that we maintain a <C l/mc 
but a ^> g 2 T/m 2 D . The latter condition is essential to make sure that "bad" hard thermal 
loops arising from modes with k ~ 1/a, which have the wrong dependence on k, uj and in 
particular are not rotationally invariant [p6fl , are subdominant to "good" hard thermal loops 
due to the particles. The other way of taking the limits is to make mo 3> so the wave 
number (1/a) falls between the gT and g 2 T scales. 

In each approach, the model correctly generates the effective HTL dynamics up to power 
corrections in (g 2 aT), (Q/g 2 aT), g 2 T/rri£), and (g 2 T) / \am 2 D ) , except perhaps for the behavior 
of the parallel transporter. The contribution to the parallel transporter from IR fields is 
correct because for a gauge field of wave number k 1/a the replacement of the straight 
line path with the "jagged" path actually used (see again Figure |3|) gives the right behavior 
up to corrections suppressed by 0{k 2 a 2 ). The question is, how does the UV contribute to 
the parallel transporter? 

In the case where we make vtld ^> 1/a, all of the lattice degrees of freedom have over- 
damped evolution. None of them propagate and there are no problems from hard thermal 
loops arising from classical lattice degrees of freedom with k ~ 1/a. Except for the parallel 
transporter, the systematics are then best under control. However the log appearing in the 
parallel transporter is obviously cut off by the inverse lattice spacing, log(mi)/ ' g 2 T) becomes 
C + log(l/ g 2 aT) with C a coefficient to be determined by matching. Since we know that the 
lattice regulator is not rotationally invariant we expect that C will be direction dependent, 
leading to a rotationally non-invariant correction. However, this correction is only at next 
to leading order in an expansion in log(l / ' g 2 aT). 

In the other case, < 1/a, all of the physics which sets both upper and lower limits 
on the log in Eq. (^) lie at k <C 1/a, and naively we should get the correct behavior at 
leading order in g. We do not, however, because the most UV lattice modes have dispersion 
relations which "turn over," k > uj. Diagram (a) from Figure [| can occur with the emitted 
gluon on shell. In fact this is just Cerenkov radiation, which is permitted because the group 
velocity of the most UV lattice modes is subluminal, while the hard particles move at the 
speed of light. Hence, while there is no contribution to the parallel transporter from modes 
with 1/a ^> k ^> mp, there is a contribution from k ~ 1/a, and again we get a rotationally 
non-invariant correction, \og{mr,/ g 2 T) becomes C + log(m£,/g 2 T). 

To determine the value of C we need to repeat the arguments leading to Eq. ( P0|) but 
using lattice gauge fields multiplied along a "jagged" path which stays closest to the straight 
line path under consideration. For a Wilson line of length I = La in a direction with unit 
vector pi, with pi in the first octant, the equivalent of Eq. fl20|) is 
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and the values in other octants follow from cubic symmetry. Here the sum over e is over 
a basis of the two unit vectors satisfying J2i e i sin(fcja/2) = 0. The integral over (x,y,z) 
appearing in the definition of F averages over starting positions for the Wilson line within a 
lattice cell, and square brackets always mean the argument is rounded down to an integer. 
The rounding down to an integer makes it very difficult to evaluate the expression analyti- 
cally. It is also tricky to find p(k,u) on the lattice, away from the mr> > 1/a or mo <C 1/a 
limits. 

In the infrared, meaning k <C 1/a, the sum in F is well approximated by an integral, 
recovering Eq. fl20p . In the UV, while Eq. (|31|) is manifestly cubic symmetric, it is not 
spherically symmetric: the k integral is over a cubic range, and neither F(k,u) nor p(k,u) 
are rotationally invariant. The integral has been computed by Arnold and Yaffe in the large 
L limit for the special case that pi lies in a lattice direction and the spectral density is 
concentrated at uo = P5 |, although they did not present their calculation in this context. 



I will now compute it in the opposite limit, <C 1/a, but still along a lattice direction. If 
L is large and Pi lies in a lattice direction (say, the z direction) then F(k,u) ~ 2nl5(k z — 
uj mod (27r/a)). For m^ < 1/a the dispersion relation is the free lattice dispersion relation 
uj 2 a 2 = Hj(2 — 2 cos(kiCi)) and the spectral density p all lies on shell, and the integral reduces 
to 

2 J-n/a (2tt) 3 (k 2 ) 2 V V a J ' K J 

where k{ = 2 sin(A; l a/2)/a. Evaluating numerically gives .0466539iVg 2 /T, which is .586 times 
the coefficient of \og{mn/ g 2 T). This demonstrates that there are extra UV contributions 
arising from k ~ 1/a in the case mo <C 1/a, which though not terribly large are not 
negligible. 

The upshot is that however mo is tuned, the UV lattice degrees of freedom will introduce 
a rotationally non-invariant contribution to the radiative correction to the Wilson line which 
determines the parallel transport of a particle's charge. This spoils the rotational invariance 
of the IR HTL effective theory at next to leading order in log(l/<?). The best we could do 
would be to average the value of Eq. fl3lD over angles. Then an appropriate choice of lattice 
spacing a in the regime with mo 3> 1/a could match the value of the log(m£>/ g 2 T) as closely 
as possible to the quantum theory value. This is the best option I see for determining the 
sphaleron rate beyond leading order, but it does not eliminate all systematic errors even at 
leading parametric order in g. 

I have not shown that the same problem will arise for other possible lattice implementa- 
tions of HTL effects; but since such implementations must generally involve parallel trans- 
portation on the lattice I expect that the problem discussed here is general. Certainly, when 
proposing some other numerical implementation of HTL's, the burden of proof must lie on 
the side of showing that problems from rotational non-invariance of Wilson lines do not arise. 

4 Numerics 

Now that I have discussed the establishment of the effective theory, Eq. (13), I will discuss 
how to make a lattice model of that effective theory and how I compute the leading log 
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coefficient of the sphaleron rate by using that model. First, define a Langevin time r related 
to the time appearing in Eq. (|T5|) via 



Note that r has dimensions of length squared, not length. The effective infrared theory then 
has the familiar form 

[D T , A t ] a (x, r) = [Dj, F 3l }\x, r) + tf(x } r) , (tf(x, T )$(y, r')) = 2T6(x - y)5(r - r') , 

(35) 

which is a Langevin equation. Studying it, I will find the diffusion constant for Chern-Simons 
number 

r r ((iV cs (r)-iV cs (0)) 2 ) 

r La „ gC vin = Jim Jim — , (36) 

numerically. To do that I find the diffusion constant per lattice site per a 2 of Langevin time, 
and multiply by a -5 . On dimensional grounds r Langevin must be of order a^/T 5 . Using Eq. 
(P^), the relation between r Langevin and k' defined in Eq. (g) is 

3iv r L 

Before getting further it is worth commenting that the IR behavior of Eq. ( |35| ) is insensi- 
tive to the UV regulation and the limit in which that regulation is removed, exists. Though 
this statement appears banal, it makes Bodeker's effective theory completely different from 
the classical Hamiltonian dynamics. The essential difference is that there are no long time 
scale correlations for the UV fields in the Langevin evolution; a mode with wave number 
k gets randomized in Langevin time r ~ 1/fc 2 , which is much faster than the natural time 
scale for the evolution of IR fields. Hence the IR fields see the average over all excitations 
of the UV fields. The influence of the UV must be purely thermodynamic, and we know 
from the super-renormalizability of 3-D Yang-Mills theory that the thermodynamic influ- 
ence of the UV on the transverse sector is well behaved. For a more rigorous presentation 



of the argument see |T{|, who show that any purely dissipative update algorithm will give 
a good continuum limit. For the Hamiltonian system, on the other hand, UV modes are 
propagating; the unequal time correlator behaves like cos(kt) rather than like exp(— k 2 r). 
The existence of long time scale correlations of the UV fields is what makes the HTL effects 
important to the IR dynamics. Because of this difference, we can expect a good small lattice 
spacing a limit to exist for the Langevin time dynamics, and it is worth it to try to control 
systematic errors. 

I discuss the continuum r, spatial lattice implementation of Eq. (|35D in the appendix; 
here I will just mention how I discretize the time update. I define a time step At = a 2 A, A 
a pure number much less than 1. The fields will be well defined at times nAr, n an integer. 
The noise is constant in each interval [nAr, (n+l)Ar); its value at each point, direction, and 



Lie algebra direction is drawn from the Gaussian distribution with mean value y2T/a 3 Ar; 
its value at each point, direction, and Lie algebra direction, and in each time interval, is 
independent. To determine the fields at time (n + l)Ar from the fields at time nAr I use 
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10.26 ±0.79 



Table 1: Results for n' at two lattice spacings and two lattice volumes. The results show 
excellent spacing and volume independence. 



the following second order algorithm: I compute dA/dr(nAr + 0) and use it to predict 
A((n + 1)At); then I average the values of dA/dr at the starting point and at the predicted 
end point, and use this average to update A(nAr) to time (n + l)Ar. The step size errors 
are 0(A 2 ). In particular, if I were studying the free theory a mode with wave number k 
would be updated with step size errors ~ (fc 2 Ar/2) 2 . The direct errors in the update of the 
IR fields are tiny. However, the most UV modes are only updated correctly if A <C 1. The 
UV fields influence the IR modes radiatively, so we do need A to be small; but the radiative 
corrections are suppressed by 0(g 2 aT), and it is not too difficult to get A small enough to 
make step size errors subdominant to statistics. 

Unfortunately the above update is quite inefficient. However there is a much more effi- 
cient algorithm for dissipatively updating the fields, the heat bath algorithm. Rather than 
applying a very small step of Langevin update to each lattice link in parallel, the idea is to go 
through the links of the lattice at random, performing a complete heat bath update of each 
link. The relation between Langevin time and the number of links updated is discussed in 
Appendix [A]; in particular it is possible to make a very accurate match between the Langevin 
and heat bath time scales by measuring the autocorrelations of some IR observable. I apply 
the small A limit as part of the matching, so that is taken care of. 

Appendix [A], together with previous work |29]| , shows how to control lattice spacing 
systematics so they first appear at 0(a 2 ). Also, it is possible to define a lattice measurable 
to use for iVcs which is topological and will eliminate systematic errors in T due to the 
definition of Nq§. in fact, two fairly efficient techniques are available [[11], |3(J; here I will use 
the method developed in It remains to take the large volume and time limits. It was 
shown in [|K| that finite volume systematics are negligible on cubic toroidal lattices larger 



than S/g 2 T on a side. To be doubly sure, I have used a lattice 16 /g 2 T on a side; as a check 
I measure T also on a lattice of half this size to check that the result is the same. Taking 
the infinite time limit is tied up with the problem of converting a Langevin time series for 



iVcs into a measurement of F. I use the same analysis techniques as [11 



To verify good control of lattice spacing systematics I have made measurements of T 
at three lattice spacings, a = 2/3g 2 T ((3 L = 6), a = 2/hg 2 T ((3 = 10), and a = 2/7g 2 T 
[fi L = 14). The results are presented in Table |l|. Finite volume and spacing systematics 
are under control. In particular, varying the lattice spacing by over a factor of two leads to 
corrections smaller than the statistical errors. This makes large a extrapolation unnecessary, 
which is very important, since numerical cost rises as (1/a) 5 . 

If I had used naive rather than radiatively corrected relations between lattice and contin- 
uum parameters, then the value of k! at the three lattice spacings would be 17.7, 15.0, and 
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13.3 respectively. Such strong dependence is because converting T from lattice to continuum 
units involves the 5'th power of a, and the radiative corrections are at leading order a shift 
in the meaning of a from the naive value, by roughly 10%, 6%, and 4% for the three lattice 
spaces used. At two loops there is a further shift, estimated to be of order (and probably 
less than) the square of the first order shift p9fl ; around 1%, .36%, and .2% respectively. The 
latter two are negligible compared to the statistical errors, even after taking account of the 
5'th power dependence. It is less clear how to estimate the importance of 0(a 2 ) nonrenor- 
malizable operators; but if I estimate all 0(a 2 ) errors by using the three lattice spacings to 
extrapolate to a = 0, assuming errors proportional to a 2 , the result lies within the error bars 
of the two finer lattice data, and the error in the extrapolation is dominated by the error in 
the finer lattice data. Hence I will adopt the middle lattice spacing result and its statistical 
error bars as the best estimate. 



5 Interpretation and systematic corrections 

The numerical result is that 

T = (10.8 ± 0.7) (log ^ + 0(1) ) ^-paLT 4 + (higher order in g) , (38) 
V 9 T J m D 

but it remains to determine or estimate the rate at the realistic standard model value of 
mo, m 2 D = (ll/6)g 2 T 2 with g 2 ~ 0.4. For this value, log(rriD/g 2 T) ~ 1.5, and the 0(1) 
correction may be quite important. 

In fact we might expect that the log needs to be quite large before it dominates the 
0(1) "correction". The leading log behavior is based on the hard particles propagating 
only a short distance before undergoing a collision which randomizes their charge. This 
"short" distance is 4ir/(Ng 2 T\og(m£)/g 2 T)), with N = 2 since we are in SU(2) theory. The 
leading log contribution to the free path is ~ (6/g 2 T)/ log(mx>/ g 2 T). For comparison, above 
I confirm Ambj0rn and Krasnitz' result [ffUl that a lattice only 8/g 2 T across is already large 
enough to give continuum like behavior for the sphaleron rate. The nonperturbative length 
scale characterizing baryon number violating processes must be shorter than this, perhaps 
by a factor of two. Hence the log will need to be quite big before the "short" distance really 
is short compared to the scale which is setting the physics. This supports the expectation 
that there will be large corrections to the leading log. 



5.1 estimate using Laine and Philipsen's results 

As I noted already in subsection [2.3| , the log arises from the behavior of a Wilson line, and 
the same Wilson line appears in the definition of the Debye mass beyond leading order. In 
this context the value for the 0(1) correction to the leading log behavior has been found 
by Laine and Philipsen. In that case, [log(m£>/ g 2 T) + 0(1)] has 0(1) = 6.7. If the same 
number held for the sphaleron rate, then using the standard model value of to evaluate 
the log, the leading behavior would be T = 89(g 2 T 2 /m 2 D )a\ v T A = 48a^/T 4 . This is a crude 
way of estimating the nonleading corrections, though, and I do not take it too seriously. 
In particular there is evidence that the length scale relevant for baryon number violating 
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processes is longer than the ~ ^-/g 2 T Laine and Philipsen find for the 0(g 2 T) correction to 
Debye mass; the baryon number violation rate on a cubic toroidal lattice 3/g 2 T across is 
over 1000 times slower than for a large volume (|30f| , so physics must be going on involving 
lengths at least half as long as 3/g 2 T. 



5.2 estimate using Hu, Miiller, and Moore's results 

Another way of trying to determine the subleading corrections is to use the value of the 
leading log coefficient to correct my results with Hu and Miiller. There, we used a technique 
which included hard thermal loops, but in a way which does not correctly reproduce the 
subtleties of the Wilson line responsible for the logarithmic dependence of T on mp. As 
discussed in subsection |Q| the log arises because excitations more UV than the g 2 T scale 
make the Wilson line effectively randomize the charge of a propagating particle, and there 
is a log in the reciprocal length for randomization. The best approach would be to compute 
the angle averaged value of the reciprocal length, in the quantum theory and in the lattice 
theory actually studied. Unfortunately, so far I can only compute the reciprocal length in 
the lattice theory in one direction, and only in the cases <C 1/d or ^> 1/a. The data 
in [13| are taken at m D ~ 1/a, not ^> 1/a, so it is at least reasonable do the match using 
the small m^a approximation. I will make do with the log evaluated in the one direction 
where I can do the integral, which gives a difference in logs between the two theories of 
log (mo (lattice) jrau (continuum) + 0.59. 

I will use the datapoint from that paper taken with the largest value of mp, because the 
"wrong" lattice mode induced HTL's really are strongly subdominant to the "right" particle 
induced HTL's for this case. Redoing the match between the lattice and continuum time 
scales, which was performed wrongly there because we did not have the 0(a) calculation 
performed in Appendix [A] of this paper, revises the result from T = 53 ± 5(g 2 T 2 / 'm^a^yT 4 
down to 50 ± 5 (g 2 T 2 /m 2 D )a\ v T 4: . This result was obtained at m D ~ 4g 2 T, so the difference of 
logs between the lattice theory where this number was computed and the quantum theory at 
the physical value of m 2 D is about .63 due to the values of mp plus .59 due to UV contributions 
present on the lattice but not in the continuum theory. Using the determined coefficient of the 
leading log term, I should correct the diffusion constant we found downwards by 1.22 x 10.8, 
giving 37 ± 5, with only the statistical error bar shown. (Substituting in m 2 D = llg 2 T 2 /6 
gives T = 20 ± 3a 5 w T\) 

Using k! to correct the old data in this way assumes that 

= k'o^T 4 • (39) 

We really only know that this true in the large (mj^/g 2 T) limit, where the leading log 
expansion is valid. It may have quite nonnegligible corrections at realistic values of ttid, which 
would appear in a systematic expansion in \og(rri£>/ g 2 T) as inverse powers of the log; the 0(1) 
correction to the leading log behavior would really be 0(1) = C\ + 02(log(r«D/ i Q 2 T))~ 1 + . . .. 
Intuitively I expect that the real mo dependence will be weaker than the leading log suggests 
(meaning O2 > 0), though I cannot give a cogent argument to show this is so. In this case I 
have performed an overcorrection, and the real rate would be higher. For now I will accept 
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the corrected answer as the best current guess, but I take a systematic error bar of order 
30% to cover both the rotational non-invariance not handled correctly in the correction, and 
the fact that the correction may be an overestimate at realistic mp. 



5.3 corrections which are formally parametrically suppressed 

There is a further cause of systematic error in the determination of the sphaleron rate, arising 
from corrections to the g <C 1 approximation. I will mention the two such corrections which 
I think are the largest; fortunately they have opposite sign and the optimistic can hope that 
they largely cancel. 

One problem is that the parametric argument for the importance of hard thermal loops 
assumes ^> g 2 T, and it is not clear realistically whether this is obtained [14]. There is 
evidence that it is not. In particular, while the sphaleron rate in classical Yang-Mills theory 
depends on the lattice spacing in a way consistent with the Arnold- Son- Yaffe prediction 



T oc m^ 2 oc a, the corrections to the linear dependence are large[ll]. Of course, some 
dependence is expected, since we know now that the scaling behavior should be not r oc a 
but T oc a(log(l/a) + 0(1)); but we can determine the coefficient of the log(l/a) term by 



using the results of this paper. According to Arnold |2"7|] , we can relate the classical lattice 
results to continuum ones, approximately, by replacing m 2 D ~ .Q8Ag 2 T/a. In this case the a 
dependence of T in classical, lattice Yang- Mills theory should be 

.465*;' (log [(g 2 aT)-^} + 0(1)) ^c/^T 4 , (40) 



on the lattice 



with k,' the same as the one we find but the 0(1) correction different. This formula makes it 
possible to correct the data in |TT| to remove the logarithmic dependence on a, for instance 
by adjusting the data so they all correspond to g 2 aT =1/4, which is the value for the finest 
lattice used there. I have done so, and the result is plotted in Figure £|. 

I fit the data to the form T/a = C\ + O2C1 to find the corrections to the small spacing 
limit which do not arise from the log. The fit is very good, see Figure |4], but the coefficient 
O2 is quite substantial. The physical value of HTL strength, using mr> = (ll/6)g 2 T 2 and 
g 2 = 0.4, corresponds to g 2 aT /A = .037, where the correction accounts for an 11% shift 
from the small a limit, with the actual value falling below the large mp extrapolation. The 
correction reduces the sphaleron rate. Given the other systematics in play it is probably not 
fair to call this a measurement of the correction to the mo» g 2 T limit; rather I will call it 
an estimate to tell how large the systematic error is. 

Another correction which is parametrically suppressed but not necessarily very small 
arises from QCD scatterings of quarks. Quarks are responsible for almost half of m 2 D , and 
hence almost half of all HTL effects. But quarks scatter strongly. While a strong scattering 
does not disturb a quark's electroweak charge, it does change its direction of flight, whereas 
the calculation of the HTL's is made assuming particles maintain straight line trajectories. 
If the free path for strong scattering were < l/g 2 T this would make an 0(1) correction to 
the quarks' influence on IR physics. Actually the quark free path for large angle scattering 
is parametrically order (a 2 log(l/g s )T) _1 and the correction is formally 0(a 2 /aw) ~ 9 2 \ but 
numerically this might not be small. The effect of this correction is to increase the sphaleron 
rate. It might be possible to include this correction within the context of the "particles" 
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Figure 4: Small lattice spacing extrapolation of data in pure classical lattice Yang-Mills 
theory, taken from \TL\ . The data have been corrected to absorb the leading log dependence 



on m? D oc 1/a determined here. They show a substantial linear correction to the predicted 
r oc a behavior. This is evidence of nonnegligible corrections to the parametric mp 3> g 2 T 
limit. 



approach by calculating more precisely the size and frequency of quark scatterings and adding 
them to the dynamics of the particle degrees of freedom. 

5.4 including the Higgs 

I should also mention that the evaluation of the leading log coefficient made here was within 
pure SU(2) Yang-Mills theory, without a Higgs boson. In the context of baryogenesis we 
actually need to know the sphaleron rate in the presence of at least one Higgs boson, in the 
symmetric phase and roughly at the equilibrium temperature. The Higgs field's evolution 
is not overdamped, because the hard thermal loops for a scalar field are nothing but a 
mass squared correction. Hence the k ~ g 2 T modes of the Higgs fields evolve on the time 
scale l/g 2 T, which is parametrically faster than the gauge fields. Therefore, on the l/g 4 T 
time scale, the gauge fields only see the thermodynamic average over all Higgs fields in the 
fixed gauge field background^. This could be simulated by including the Higgs field in the 
Hamiltonian, and evolving it with heat bath dynamics, but updating the Higgs fields much 
more often than the gauge fields. Then we would have to extract the limit as the Higgs field 
update is made infinitely faster than the gauge field update. I have not yet attempted to 
do this, but I anticipate that for parameters which make the phase transition strongly first 
order, the effect should be a slight reduction of the leading log coefficient for the sphaleron 
rate. 

6 Dam Son pointed this out to me. 



23 



6 Conclusions 



In the formal small g limit, in which an expansion in log(l/g) 3> 1 is justified, the sphaleron 
rate in SU(2) Yang-Mills theory is 

T = (10.8 ± 0.7) flog ^ + 0(1)) ^a^-T 4 . (41) 

The value 10.8 ± 0.7 is clean; the errors are dominated by statistics, systematics are well 
under control. 

However, interpreting this result to get the sphaleron rate at the realistic values of m 2 D = 
(ll/6)g 2 T 2 and g 2 ~ 0.4 is very problematic, because the 0(1) correction is not subdominant. 
The reason is that the expansion in large log(l/g) corresponds to treating the length scale 
2n/[g 2 T\og{mo/g 2 T)] as much shorter than the scale relevant for baryon number violating 
processes, which is around ~ 4/g 2 T. The best current estimate for the 0(1) correction comes 
from using the result of this paper to correct previous results of Hu, Miiller, and myself, which 
effectively were using the wrong value for the log. The result at the physical values for mo 
and g 2 is F ~ (20±3)af l/ T 4 , with only statistical errors quoted. Unfortunately the correction 
procedure is not well under control; I estimate that the systematic errors should be taken 
to be at least twice as large as the statistical ones. There are also corrections to the g <C 1 
approximation, probably at around the 10% level. 

Besides making it possible to correct previous results, the leading log approximation and 
Bodeker's effective theory provide a clean test-bed for determining how good an approxi- 
mation it is to treat the symmetric phase of Yang-Mills Higgs theory using just Yang-Mills 
theory. Checking how much difference the Higgs field makes is an interesting project for 
the future. It would also be interesting to study the dependence of k! on the number of 
colors, since the sphaleron rate in the strong sector is also phenomenologically interesting 
for baryogenesis. 
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A O(a) lattice time scale renormalization 

In this appendix I compute at O(a) the lattice to continuum correction to the Langevin 
time scale, and the relation between the continuum Langevin time scale and the amount of 
heat bath applied. The calculation of the Langevin time scales will be quite technical and 
will depend to some extent on previous work relating the lattice and continuum theories at 
the thermodynamic level, found in J^J. I will relate the heat bath depth to the continuum 
Langevin time scale by making a nonperturbative lattice measurement comparing it to the 
lattice Langevin time scale, and then using the analytic relation between this and the con- 
tinuum Langevin time scale. It would also be possible in principle to make the connection 
directly, but the analysis is difficult and I am lazy. 
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A.l lattice and continuum Langevin time scales 



Reference [£9| explains how the lattice to continuum relation for thermodynamical proper- 
ties of 3-D gauge or gauge-Higgs theory can be studied as an expansion in g 2 aT, where a 
is the lattice spacing and g 2 is the coupling in 4-D notation. The combination g 2 T is the 
coupling constant of the 3-D theory. Since it is dimensionful, a perturbative matching be- 
tween continuum and lattice theories must be an expansion in g 2 Ta on dimensional grounds. 
In Yang-Mills Higgs theory with fundamental or adjoint scalar fields the leading terms can 
behave as 1/a because the theory contains dimension 2 operators, but in strict Yang-Mills 
theory the lowest order operator is dimension 4 and the leading corrections are 0(a). Fur- 
ther, the only 0(a) correction is a rescaling of the coupling, equivalent to a rescaling of the 
physical length scale, and it only arises from one loop diagrams, and has been computed 



Arnold, Son, and Yaffe have recently demonstrated that, because of super- renormaliza- 
bility and general considerations worked out by Zinn- Justin and Zwanziger the same 



applies to the dynamics under Langevin dynamics In particular, any high dimension 



corrections which can appear in Eq. flT5"D would change the dynamics at 0(a 2 ), simply 
because there are no gauge invariant, P even dimension 5 operators in Yang- Mills theory. 
The only possible 0(a) correction, besides the thermodynamic one already mentioned, is a 
rescaling of the Langevin time scale, which will arise exclusively at one loop. This opens 
the possibility of computing the 0(g 2 aT) corrections between the continuum and lattice 
Langevin dynamics, by performing a one loop match. It also means the Langevin dynamics 
can be replaced by heat bath dynamics at the cost of another 0(g 2 aT) shift, which I will 
measure, rather than compute, in the next subsection. 

In the continuum theory and in temporal gauge the Langevin evolution is 

The Langevin equation for the lattice theory must be written in terms of the unitary parallel 
transporter matrices Ui(x) and their derivatives. By definition the matrix Ui(x) G SU(iV) is 
the matrix such that, if $ is a fundamental representation object transforming as an object 
at point x + ai, Ui(x)<& is the parallel transport to point x; thus Ui(x) "lives on" the link 
between the site x and the site x + ai. Writing D\ for the left acting derivative, 

D a L U = igaT a U, D a L F(U) = F(U -> U + D a L U) - F(U) , (43) 

the Langevin equation for U is 

^ = -Dia t (x)(p L TDlH K , + £(*)) , (&{z)${y)) = %§#^(r " r') . (44) 
aT g a Pl 

Here T a is a fundamental representation Lie algebra generator with the standard normal- 
ization; i^KS is the Kogut Susskind Hamiltonian, the sum over elementary plaquettes of the 
trace of the plaquette, 

H KS = J2±-^U n ; (45) 
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and (3l is the inverse temperature in lattice units. At leading perturbative order (3l = 
A/(g 2 aT), but it receives radiative corrections, computed in ||29|| , which shift it by a constant. 
What I write here as (3l is /3 na i ve in the notation of [p9| , but in the body of the paper I 



have always used the 0(a) improved definitions when I refer to a or /3l- The combination 
PlHks equals H/T of the continuum theory, up to radiative corrections and high dimension 
operators which correct infrared physics at 0(a 2 ). The radiative corrections are absorbed 
up to errors which are 0(a 2 ) by the one loop radiative correction to /3l, which I will present 
in due course. 

Our task is to compute to 0(a) the relation between the Langevin time scale r in the 
continuum and lattice cases. The correction at loop order I will be 0(g 2 aT) 1 , with the g 2 T 
from loop counting and the a to balance dimensions; so we need only go to one loop order. 
Part of the correction is from the shift in {3 already mentioned, but there are also corrections 
from the relation between U and A and from radiative differences between the lattice D a and 
the continuum derivative. I will be satisfied to perform the calculation in a particular gauge, 
strict Coulomb gauge, without checking for the gauge independence of the result. By strict 
Coulomb gauge I mean that at every Langevin time the 3-D configuration is in 3-D Landau 
gauge. I should also fix a global time dependent gauge ambiguity; but this is irrelevant at 
the level of perturbation theory. 

In Coulomb gauge the U matrices are all close to the identity matrix, with a departure 
of order f3 L ; so it makes sense to define a lattice gauge field A\ through 

Ui(x) = exp(igaT a Af) . (46) 

the A field so defined has the same normalization as the continuum A field at tree level, but 
there are radiative corrections. The value of the radiative correction to A is obtained in the 
theory with an added scalar field by matching the one loop values of the gauge-scalar three 
point vertex at small transfer momentum. There is a contribution from one loop vertex 
corrections and from one loop scalar wave function corrections, and neither depends on the 
number of scalar particles in the theory, so the answer is the same in the pure gauge theory. 
The one loop correction was computed as a byproduct in |2D| in a general gauge, for SU(2). 
The result, in Coulomb gauge but for SU(iV) gauge theory, is 



A" (continuum) 



\ 4 \ 18 4tt 4vr , 



A a { (lattice) . (47) 



The numerical constants £ and £ were first defined in ||32|| , and their numerical values are 
£ = 3.175911536 and £ = 0.152859325. 

Next I need to find the relation between applying the lattice and continuum derivatives. 
For this purpose it is actually more convenient to write the update rule in terms of center 
acting derivatives, 

D a c U = U 1/2 (igaT a )U 1/2 , (48) 

because the formulation will then be parity symmetric. Naively one would expect that if we 
set dUi/dr = E^D^U, that dAf/dr = E®, and indeed this is correct at leading order in A. 
But beyond leading order there are corrections; expanding both sides gives 

dUi d . OT1 „, dAI g 2 a 2 T b T c ( Ah dAS „ A 
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ig a 



3 „3rpdrperpf 



E*D a U 



1 + igaT a A° 



( dA f \ 
[AfA?-^- + permutations J + 0(a 4 ) 

2 2 

9_^ T b T C A b A C + 



2„2 



(49) 



igaT d E? x 



1 + igaT e A\ - ^T f T h A f i A'l + 0(a 3 ) 



(50) 



We must expand to 0(a 2 ) corrections because (A 2 ) ~ 1/a. Equating the expressions, after 
some work we obtain 



dA^ . , . 
— — (lattice) 
dr 



E" 



hid 
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f abe f cde AlA'r + 0(a 3 



(51) 



To save some writing I left off marking that all A, E above are in position space and at the 
same coordinate position. The above correction is strictly a lattice effect and the equivalent 
continuum relation is dA^/dr = Ef. 

Naively I should now equate Ef with —(3lD^;Hks + ^, evaluate the (UV dominated) mean 
value of the correction between the lattice and continuum relations, and thereby determine 
the rescaling of the Langevin time. This would not be strictly correct, though, since the 
Langevin equation, Eq. (|4]), gives time evolution in temporal gauge. The temporal gauge 
evolution breaks the Coulomb gauge condition, and if we are to do the calculation in Coulomb 
gauge we must make a time dependent gauge change to maintain the Coulomb condition 
at all Langevin times. There is some danger that the gauge changing will also lead to a 
radiative correction to the time relation; it turns out it does not, but I will go through the 
calculation anyway to show that it does not, since the cancellations may be special to strict 
Coulomb gauge. 

The Coulomb gauge condition is Dl ■ A a (x) = 0, which means 



£A a (*)-A a ( 



X 



at 



0. 



(52) 



This does not look cubic invariant because our labeling associates U^, and hence A iy on the 
link between x and x + ai with the site x; it might better be associated with the point 
x + cm/2, which would make the cubic invariance more obvious. 

To maintain this Coulomb condition, we must apply a time dependent gauge trans- 
formation at each site. The difference between Coulomb gauge and temporal gauge sat- 
isfying the Coulomb gauge condition at r = will be a gauge transformation by A = 
Texp / igT a G a (r)dT), where the T means the exponential should be time ordered with re- 
spect to t. The value of G is fixed by the requirement that the Coulomb condition remain 
true; we must choose G so the departure from the gauge condition due to evolution of the 
fields and due to the action of G cancel. In an infinitesimal time interval cir the gauge change 
alters U through 



SUi(x) = (1 - igT a G a (x)5T)Ui(x)(l + igT b G b (x + ai)5r) - U^x) 
Expanding U on each side, the G contribution to the time evolution of A is 



(53) 



dA a { 

~dT 



(from G) 



G d (x + at) - G d (x) 



5, 



ad 



12 



fabefcdeA\{x)A\{x) - 
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G c (x + ai) + G c (x) 



-gf ahc A\{x)^ 1 ^ . (54) 

Now we need to substitute this expression, and the relation between E and dA/dr, into 
Eq. ([52]) to determine the relation between E and G. The result, Fourier transformed to 
momentum space, is 

• 2 2 

= ikEf(k)- tJ ^kfabefcdef ^{l)A^(m)Ef{k - I - m) - 

24 Jim 



-k 2 G(k) - igf abc k l jcos(al l /2)A h l (k - l)G c (l) - 



-^rkkfabefcde [ A^A^k - I - m)G d {m) . (55) 

IZ Jim 

Here summation over vector indices is implied in terms where the index appears at least 
2 times. I have used the conventional shorthand ki = (2/a) sin(fcja/2), k 2 = J2ik 2 , and 
Ji = J d 3 l/(2n) 3 , where the range of the integration is [— n/a, n/a} 3 . Indices are summed 
whenever the index appears an even number of times, which can be more than twice in 
cubic symmetric but rotationally nonsymmetric expressions; but if an index appears an even 
number of times on each side of a + or — sign I am re- using it, it is an independent index 
in each expression. To get the continuum version of this expression, drop the 0(a 2 ) terms 
and set k = k, cos(/ca/2) = 1. 

We can determine G perturbatively in g by expanding Eq. (|55"D in powers of g. The 
result to 0(g 2 ) is 

Le? „ rkd 



G a (k) = i^ + gU c j^ico S (ka/2)A^k-l)E^l)- 



k 2 Jim 

-ig 2 f ab ef cd e j cos(W2) cos(m ia /2) A\{k - l)A%l - m)E d k {m) 

Jim k 2 l 2 m 2 

ig 2 a 2 ki 



2 2; 

""" '\U.U I Al(l)A^m)Ef(k - I - m ) 

J lm 
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W 2 cl 2 f kikifh 



f ab ef cd e / -±^A$(l)A${k - I - m)Ef(m) . (56) 



12 Jim k 2 m 2 

This in turn must be substituted into Eq. (|54j ) to find the true value of dA/dr in Coulomb 
gauge, 

dA ^ = U.^ME^igf^J^ 



+g 2 2 fabefcde 

J lm 



m) — 



1 / mim^bij kikjkjrrik 



A){k - I - m)A c j (l)E d k (m) - 
5ij ) cos(/ja/2) cos(mfca/2) x 

(57) 



12 \ m 2 k 2 m 2 

IhJft'l I r k{kj 



l 2 m 2 \ k 2 

A b Jk - l)A c k (l - m)Ef(m 
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Below I will be interested in the case where E is uncorrected with A, in which case the 
mean value of the term linear in A vanishes and we can substitute the leading order Landau 
gauge A field correlator, 



(A1(k)A b Al)) = 6 ab 6(k + l)[8i 



(5? 



into the remaining terms. The contributions from the terms with (1/12) and (1) in front 
vanish, only the term with the leading (1/24) coefficient contributes. It requires that we 
perform the integral 

a 2 [L = a— , (59) 
Ji P 4vr ( j 

which is the definition of the constant S which appeared earlier. We also need the integral 



P 



2\2 



The result is 



dA-{k) 
dr 



(lattice) 
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9 4 4vr 



(60) 



Em 



(61) 



valid when E is uncorrelated with gauge fields. The first factor projects out the trans- 
verse component of E and is responsible for maintaining Coulomb gauge. The continuum 
expression is the same but with a set to zero. 

Using the previously established relation between lattice and continuum gauge field nor- 
malization, the relation for the continuum normalized gauge field is 



cL4^(£;)(contin) 
dr(latt) 



4 V 6 47T 47T 



(62) 



valid for the gauge field response to the random force part of the Langevin equation. I will 
not attempt to study the response of A from the Hamiltonian gradient part of the Langevin 
equation, since it would involve understanding the radiative corrections to D a H and would 
lead to an E which might be correlated with A fields. Since we know that the Langevin 
equation correctly thermalizes the lattice system when we use the radiatively corrected value 
of (3l, it is sufficient to study the response to the random force alone to determine the rescaling 
of the Langevin time scale. 

Note that the correction appearing in Eq. (|62|) is precisely one quarter of the radiative 
wave function correction for an adjoint scalar field in 3-D lattice gauge theory when the scalar 
self-coupling vanishes, see ||29|| . This is perhaps not too surprising. In the real time theory, 
the time evolution is generated by the electric fields, which appear in the thermodynamics 



as the Aq field, an adjoint scalar with zero self-coupling. In [|ll|] Turok and I speculated 
incorrectly that the time scale correction would not contain any large tadpole corrections; 
but Eq. (|62"D is one quarter the adjoint scalar wave function correction, which does contain 
tadpoles. Our incorrect guess was based on analyzing the abelian theory, where the Aq field 
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turns out to be a free field. The correction found here vanishes in the abelian case, for 
compact or noncompact implementations. 

Now I will finish the relation between time scales. I have just shown that the response of 
the A field to the random force, shifting its normalization to correspond to the continuum 
theory normalization, is 

Hir< 1+ »^ < 63 > 

where (1 + corr) is the quantity in brackets in Eq. (0). But the autocorrelator of £ is not 
the same as it would be in the continuum, because of the radiative corrections to (3l- The 
mean square change in Af over a Langevin time Jr(latt) is, using Eq. (f4*4|), 

<(^(*)) 2 > = (1 + corr) 2 2m( 3 latt) - 4 . (64) 

However, (3l = 4/g 2 aT only at leading order in (3l- Beyond leading order, in the pure gauge 
theory, it is [^9] 

4 9 4 V 3JV 3 47r/ V ; 

The continuum theory Langevin equation would cause a mean square change to A of 



2T<5r(contin) 



mi (x)Y) = , (66) 

Lb 



so the relation between time scales is 
<5r(contin) 



5r(latt) 



g 2 aT fN 
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(1 + corr)% = 1 + 2—- — — - — -± — — . (67) 



The numerical value of this expression for N = 2 is 1 — 0.3189(g 2 aT/4). 

I have now related the lattice and continuum Langevin time scales at 0(a). It is also easy 
to show that, for the Hamiltonian system, the correction between the lattice and continuum 
time scales is exactly half as large. However this is less useful in that case because, while this 
correction is technically correct for determining the time falloff of correlators over very short 
time scales, the IR dynamics on longer time scales receive HTL corrections which depend on 
the lattice spacing as 1/a + 0(1). For the technique of Hu and Miiller, HTL's are included 



by adding "particle" degrees of freedom. In |13| we work out the correction for time scales 
in the limit m 2 D 3> g 2 T/a, in an approximation which corresponds to setting (1 + corr) = 1. 
The correction found here changes our result there from being ti att /t C ontin = Z~ 2 to being 
Z~ 2 {1 + corr) -2 . I used this correction in Section 



A. 2 lattice Langevin time and depth of heat bath 

The effect of the heat bath algorithm on the infrared degrees of freedom is equivalent to 
Langevin evolution. I will show this at tree level, which gives a tree relation between Langevin 
time and the number of heat bath updates applied. The only 0(a) correction to this relation 



30 



possible is an 0(a) shift in time scales between the algorithms; to find the magnitude of the 
shift I make a direct numerical measurement. I do not attempt a perturbative calculation 
of the relation between continuum Langevin time and amount of heat bath applied, beyond 
leading order. 

First I will show that the effect of the heat bath algorithm on the IR degrees of freedom 
is equivalent to the Langevin algorithm, and I find the relation between the Langevin time 
scale to the number of links updated by heat bath, at leading order. The way the heat bath 
algorithm works is as follows: 

1. pick a link on the lattice at random. 

2. replace the connection U on the link with the one which minimizes the Hamiltonian. 

3. multiply this link by a random SU(2) element chosen from a distribution centered on 
the identity, with a weight function dependent only on the arc length from the identity, 
not the direction. The weight is chosen to correctly reproduce the thermal ensemble 
on this link holding others fixed; it is approximately but not exactly Gaussian. 

This is a heat bath update; for a more precise description see [j33l| . The first part, the 
quench of the link, serves to perform the —dH/dA part of the Langevin update, and the 
multiplication by a random SU(2) element reproduces the noise part of the Langevin update. 

To see the relation between the heat bath and the Langevin update at leading order it is 
sufficient to consider the linearized theory, that is, to expand the Hamiltonian to quadratic 
order in the gauge field A. On an N x N x N toroidal lattice the relation between the 
connections U, the gauge field Ai(x), and the Fourier transform of the gauge field A(k, s) 
(s a polarization index) is (writing all Lorentz and group index sums explicitly, there is no 
implicit summation convention in what follows) 

Ui(x) = exp(Y / ^9aT a At(x))^l + Y / ^9aT a A^x)- 9 ^Y, A n^(x), (68) 

\ a / a "a 

A a (k,s) = N- 3/2 J2^i(s,k)A^(x)exp (-ik ■ (x + ai/2)) , (69) 

x,i 

A$(x) = (longitudinal piece) + N' 3/2 ^ A a (k, s)e;(s, k) exp (ik ■ (x + cm/ 2 )) > ( 70 ) 

k,s 

where £j(s, k) is a transverse polarization vector, satisfying 

^2 ei(s, k)ei(s', k) — 5 s y , (71) 

i 

and 

= 0; (72) 

i 

there are two such states for each k. The sum over k includes all k of form (2n/aN)n, with 
n a triple of integers each in the range < n, < N. Only the transverse fields are of interest 
here, at the order we are working the longitudinal part is pure gauge. Its behavior depends 
on our choice of gauge fixing, and has no influence on the transverse fields (or on physics). 
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(This is not true at higher order in g, where the longitudinal terms are responsible for the 
Fadeev- Popov determinant.) At the level we are working the value in thermal equilibrium 
of A is given by 

((A a (k, s))*A b (k', s')) = S k>k/ 5 ab 5 SiS S (73) 

k l 

Now let us analyze how the fields change when the heat bath algorithm is applied to a 
link (x,i). The terms in the Hamiltonian containing the link x,i are 



H 1 

7^ D 



T " a?T 



£ [i (A](x) + A«(x + aj) - A](x + ai) - A«{x)f + 



A (-A«(x - aj) + Af{x - aj) + Aj(x - a(i + j)) - A^x)f 



(74) 



and "quench" part of the heat bath algorithm will replace A®(x) with the value which 
minimizes this expression, 

A?(x, after) = ^ (A^x) + A?{x + aj) - A a 3 {x + ai)- 

-A«{x - aj) + At(x - aj) + A%x - a(i + ]))) . (75) 

Using Eq. ([F(]) and Eq. (|72|), and adding a term £ a to represent the noise which will be 
added, this is 

a 2 k 2 

A?(x) (after) = A$(x, before) + C ~ E ^V2 e *( s ' k)A a {k, s) exp(ik ■ (x + ia/2)) . (76) 

k,s 

Now we should compute the size of £ a . Because the Hamiltonian is expanded only to 
quadratic order, the noise is Gaussian, of amplitude set by the size of the quadratic in 
Ai(x) term in H, which from Eq. (|T4]) is (2/a 2 T) J2 a Af(x)Af(x). The amplitude of the 
noise £ a is then (no sum on a) (£ a £ a ) = a 2 T/4. 

Next we will see what impact this update has had on the Fourier mode A(k, s). Combining 
Eqs. (|9D and fl7|), we find 

A a (k, s, after) = A a (k, a, before) (l - ^^f^ ) + exp {-ik ■ (x + ai/2)) f 



4iV 3 / N 3 / 2 

AN 3 



^ aH^s kUs'J) ^ s , before) ^ ^ _ iy{x + a y 2) j (J7) 

(s',lMs,k) 



The (k, s) term in the sum is removed and included instead in the first term. It is responsible 
for the damping term in the Langevin equation. Both the noise term £ and the final term, 
which I will call the "cross-talk" term, are responsible for the noise term in the Langevin 
equation. 

To measure the magnitude of the Langevin damping term, we must compute the corre- 
lator (A*(k, s, before) A(k, s, after)). Because {A(k, s) A(l, s')) = = (A(l, s') £), we get 

(A*(k, s, before) A(k, s, after)) = (A(k, s, before) 2 ) f 1 - ) . ( 78 ) 

y AN 3 J 
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It is also important to make sure that the mean square value of A a (k, s) is unchanged by 
the update, which is the requirement that the noise have the right amplitude. Here we get 
a little surprise; squaring Eq. (ff7|) , 

after)) 2 } = (l - ^^T^) <W*, *, before) 2 ) + ((C) 2 ) 



+^ E (A a (l,s')A a (l,s'))e%s,k)e 2 (s\l)(^) 



(79) 



Using Eq. (|73|) and taking A^ 3 ^> 1, this becomes 

((A%k, s, after)) 2 ) = ((A a (k, s, before) 2 ) + ^f*' k) i-\ + \ + ^ £ e 2 (s', Z)a 2 / 2 

(80) 

At leading order in large N, the sum gives 4A^ 3 . Therefore the mean size of A a (k,s) is 
unchanged, which means that we have the correct amount of noise. We see that fully half of 
the noise actually arises from "cross-talk" between the mode of interest and extremely UV 
modes, with the other half arising from the noise explicitly appearing in the algorithm. 

To be Langevin the noise must have zero unequal time correlation. This is the case for 
£ by explicit construction, but we need to check it for the "cross-talk" noise. For the heat 
bath algorithm to act like a Langevin algorithm, the "cross-talk" contribution to A(k, s) 
from updating the (x,i) link must be independent of that from the (y,j) link, at least on 
averaging over the choice of (y,j) (which is indeed chosen randomly in the algorithm I use). 
This is the case; from Eq. fl77p, the correlation between the "cross-talks" is 

(aH 2 m 2 \ 

E E : 

16AT6 j e i( s > k)ej(s, k)ei(s' ,l)ej(s" ,m) 



(l,s')^(k,s) (m,s")^(k,s) 



x (A a (l, s')A a (m, s")) exp (ik ■ (x - y + a(i - j)/2)) 
x exp il ■ (x + ai/2) + im ■ (y + aj/2) 



= E TfiXr6 e ^' S ' fc ) e ^ s ' k ) e i( s 'i e j( s '' ex P [}i. k - I) ■ (x - y + a(i - j)/ 2 )) , (81) 

which suffers from a rapidly oscillating phase. The expression is typically smaller in magni- 
tude by N~ 3 / 2 compared to the corresponding term in Eq. ( |T9"D and its average over (y,j) 
is strictly zero. Hence there is no unequal time correlation in the "cross-talk" part of the 
noise. Note also that the cross-talk is very strongly UV dominated, which means that there 
will be no hidden correlations in the IR effective evolution because of it, at least at 0(a) and 
probably higher; it is also fortunate because the UV is most quickly randomized. 

What we have shown is that applying the heat bath update is equivalent to damping 
the A fields and applying noise. In particular, its influence on the IR degrees of freedom is 
equivalent to that of Langevin dynamics. Applying the heat bath algorithm to many links in 
succession, the rate at which a mode is damped (and the amount of noise it receives) is given 
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by Eq. (fZSf) (and Eq. (|80D ) after averaging over the direction i (since each direction is bathed 
with equal frequency). Using J2i£ 2 ( s , k) = 1, we find that it takes 12iV 3 heat bath updates 
to perform the equivalent of a 2 of Langevin update. This gives the tree relation between the 
algorithms. We can then define a heat bath time in terms of the number of links we have 
updated, r(heatbath) = a 2 (number of links updated) /(12iV 3 ). At leading order in a weak 
field expansion this is the same as Langevin time but we expect subleading corrections. 

Now we must push the analysis beyond tree level. Since the influence of the heat bath 
algorithm on the IR degrees of freedom (in fact, all degrees of freedom) "looks like" Langevin 
dynamics, the analysis of applies; up to high dimension operator corrections, which by 
power counting appear first at 0(a 2 ), the algorithms are related, in the presence of interac- 
tions, by a rescaling of all parameters. Since each algorithm gives correct thermodynamic 
behavior (after the 0(a) correction already discussed is applied), the only remaining correc- 
tion would be a rescaling of the time scale, which must be at worst 0(a) since it vanishes as 
g 2 — > (in which limit the calculation just presented is exact), and on dimensional grounds 
any 0(g 2 ) correction must be 0(g 2 aT). 

We could in principle determine this 0(a) correction by an analytic computation, ex- 
tending the one just presented to second order in g. Instead I compute the subleading 
effects by the following strategy. I choose some infrared measurable O, and measure it 
at each lattice point at a closely spaced series of Langevin times. I do the same using 
the heat bath algorithm. Then I compare the unequal time correlator or autocorrelator 
C(r — t') = (0(x,r)0(x,r')) , where the average is over the ensemble of Langevin trajec- 
tories, or in practice over coordinates and times in a single very long Langevin trajectory. 
To match the time scales, we see what rescaling of the heat bath time scale is needed to 
make the autocorrelator match the autocorrelator for the Langevin case. Any IR measurable 
will do because we know that the algorithms both behave as Langevin algorithms on the IR 
degrees of freedom, so the only 0(a) difference would be a time rescaling which will be of 
the same amplitude for any unequal time observable. 

I should explain that the reason this is worth doing at all is that, first, the heat bath 
algorithm is much faster and does not suffer from step size errors like the Langevin algorithm, 
and second, there are infrared measurables other than the topological density, for which the 
autocorrelation statistics improve much more quickly. If the latter were not true we would 
spend as much computer time making the match between techniques as it would take to do 
the measurement of T by the Langevin method. 

The measurable I choose is a fundamental representation Wilson loop after some amount 
of cooling, specifically a 4 x 4 square Wilson loop after r = 3.125a 2 of gradient flow cooling 
under i^KS- This is an infrared measurable because such a large Wilson loop samples mostly 
the infrared gauge fields, and because the cooling removes most of the UV fluctuations 
anyway. 

Incidentally, it is not too hard to compute the leading order perturbative prediction for 
this quantity. The mean trace of an Z x Z Wilson loop in SU(iV) gauge theory after a length 
r of gradient flow cooling, neglecting lattice artifacts, is 

AT Trr fAT 2 U 9 2 T f d 3 k 16sm 2 (kJ/2)sm 2 (k y l/2)k 2 x + k 2 y 

N - TrUixi = (N - 1)— — I 77T -^ 2 — y - exp(-2k r) + 0(g ) , 

(82) 



4 J (2vr)3 klk 2 y 
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In the t — > limit the integral has logarithmic UV divergences but for finite r it has a well 
defined value and is dominated by the infrared, k < r -1 / 2 . Of course, for the Wilson loops 
under consideration here, perturbation theory will be unreliable because the length scales 
involved are close to the scale where perturbation theory breaks down completely. It might 
be interesting to see whether the infrared fields are stronger or weaker than at leading order 
in perturbation theory, though. 

I measured the same site, unequal time Wilson loop correlator by measuring each 4x4 
Wilson loop in an even sublattice every a 2 / 2 of Langevin time, for a series of Langevin tra- 
jectories each about 200a 2 long, with 50a 2 Langevin time between trajectories to eliminate 
correlations between trajectories. For each trajectory I determined the autocorrelator aver- 
aged over volume and time. The autocorrelation function C(Ar) looks something like an 
exponential tail but is not well fit by one; there is some small amount of much longer time 
scale correlation caused by the slow evolution of the most infrared gauge fields. To compare 
the Langevin and heat bath time scales I averaged C(Ar) over data sets for each update 
method and determined the rescaling of the Langevin time which minimized the difference 
between the results, 

X = (C Langcvin (Ar) -rcChoat^Ar)) 2 ^, (83) 

Jo 

where I allow a rescaling both of the autocorrelation time and of the overall magnitude of 
C . I chose Ar max to be enough that the autocorrelator has fallen about 1 — e _1 of the way 
to its large Ar limit, but the result is quite insensitive to the specific choice. The coefficient 
r T at the extremum of x 2 gives the rescaling of the Langevin time scale. The multiplicative 
rescaling of C is necessary because there are small, very long time scale correlations in the 
measurable which can effectively shift one data set somewhat with respect to the other. As 
a check I compared the first half of the Langevin data I took with the second half. I find 
that r T = 1 within a small tolerance the same size as the jackknife error bars, but rc differs 
from 1 by a few percent. 

For (4/g 2 aT) =61 ran each update procedure on a 24 3 lattice for r ~ 9000a 2 . The 
rescaling of the time scales was 

Ar heat = 1.098 ± .007Ar Langevin , (84) 

with the error bar determined by the jackknife method. 

The Langevin step size used here was A = .05. The definition of A and the second order 
algorithm I used are in Section |4| in the body of the paper. To check for step size errors I 
evolved a trajectory for half as much Langevin time, using A = .025. The rescaling between 
this trajectory and the heat bath was 1.099 ± .010 and the rescaling between Langevin 
evolutions with the two step sizes was 1.001 ±.009, so Langevin step size errors are negligible 
at A = 0.05. 

I also measured Nqs during the A = 0.05 Langevin trajectory; the ratio of the diffusion 
constants for the heat bath and Langevin algorithms was 

heat (unrescaled) = 1.114 ± .058 , (85) 

J- Langevin 
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which is compatible with the measured difference in time scales, but with much larger error 
bars. The autocorrelator of the Wilson loop develops good statistics more quickly than the 
diffusion constant for Acs- 

Incidentally the mean value of the Wilson loop trace was (2 — Trf/ 4X 4) = .261, while the 
leading order perturbative prediction is .096. The infrared of Yang-Mills theory has more 
excitation than leading order perturbation theory predicts, by quite a bit on the scale of a 
square Wilson loop 8/(3g 2 T) on a side. 

For (4/g 2 aT) = 10 I used a 40 3 lattice but only r = 2000a 2 . The rescaling of the time 
scales was smaller as expected, 1.065 ± .016. Assuming the subleading correction to be 
purely O(a), we would have guessed from the (3l = 6 result that the rescaling factor would 
be 1.059 ± .004, which is within error. The rescaling at both lattice spacings are within error 
of being TLangevm/Theat = Z g , and I speculate that this is the correct analytic relation at 0(a). 
For the finest lattice spacing, I have simply extrapolated the medium spacing data assuming 
a pure O(a) form for the correction. 

The deviation of the mean value of the Wilson loop measurable from perturbation theory 
is also smaller on the finer lattice; the value is .110, while the perturbative estimate is .058. 
This is also expected, since a 4 x 4 Wilson loop is smaller in physical units in this case, and 
perturbation theory works better as the length scale becomes smaller. 
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